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ABSTRACT 

Context. The formation of the first stars (Population III; PopIII) marks the end of the dark ages of the universe, a subject of lively sci- 
entific debate. Not (yet) accessible to direct observations, this early stage of the universe is mostly studied via theoretical calculations 
and numerical simulations. An indirect window is provided by integrated present day observables such as the metal abundance or the 
diffuse extragalactic photon fields. 

Aims. We aim to derive constraints on the properties of the PopIII and low metallicity Population II (LM PopII) stars utilizing limits 
on the density of the extragalactic background light (EBL), recently derived from very-high-energy (E>100GeV; VHE) observations. 
Methods. A model calculation for the evolving EBL density produced by PopIII/LM PopII stars is presented. The model utilizes 
stellar population spectra (SPS) for zero and low metallicity stars and accounts for the changing emission of an aging stellar popula- 
tion. Emission from the dense HII regions surrounding the stars (nebula) is included. The resulting EBL density for different scenarios 
(metallicity, star formation rate, initial mass function) is compared to the limit on the EBL density. The potential for detecting a cut-off 
in HE/VHE spectra is discussed. 

Results. Assuming a maximum contribution from PopIII/LM PopII stars to the EBL density of 5 nW nr 2 s _I at 2fim a limit on the 
star formation rate (SFR) of the first stars of 0.3 to 3 M Q Mpc~ 3 yr~' in the redshift range 7 - 14 is derived. The limit depends on the 
assumed shape of the SFR and metallicity. 

Conclusions. The EBL can be used as a probe to investigate the properties of PopIII/LM PopII stars. Limits on the EBL density 
derived from VHE observations can provide constraints on the parameters of the these stars, in particular the star formation rate. 

Key words, early universe - diffuse radiation - Gamma rays: observations 



1. Introduction 

The end of the dark ages of the universe - the epoch of reion- 
ization - is a field of great interest (e.g. Barkan a~& Loeb| 2001 



|Ciardi & Ferra ra 2005 ). This epoch is associated w ith the for 



mation of the first stars (Population III; PopIII J] (e.g. Bromm & 



Larson 2004 ; Glover 2005 ), which are believed to start the reion- 
ization of the universe at redshift of about z = 10-30. PopIII 
stars form in a pristine environment, in clouds of hydrogen and 
helium with little or no heavy elements (primordial composi- 
tion). Due to the absence of heavy elements, the cooling of such 
collapsing gas clouds is likely dominated by H2 cooling through 
molecular emission lines. Numerical simulations of collapsing 
clouds with primordial composition predict very massive stars 
( 1 00-1 000 Mq) with high effec tive temperatures (~ 10 5 K) and 
short lifetimes (~ 10 6 yrs) (e.g. |Bromm et al.|l999||2002||Abel 



et al.|2002l >. Such hot massive stars produce copious amount of 



ionizing photons ( |Schaerer|2002| l and can therefore reionize the 
universe. The formation of lower mass stars is also possible, if 
e.g. the cloud cooling is driven by hydrogen-deuterium (HD) and 
atomic hydrogen (H) cooling ( Uehara & Inutsuka 2000 ; Johnson] 
& Bromm|2006 1. Other processes including turbulent fragmen- 



tation ( Klessen et al.|2005 l, magnetically-regulated fragmenta- 



tion ( Silk & Langer 2006 ) and dust cooling at very high densities 
( |Omukai et al.|2005) > could also explain stars with lower masses 
< 100 McQ 

PopIII stars produce the first heavier elements, paving the 
way for the second generation of stars. When the star form- 
ing cloud reaches a critical metallicity (Zqr ~ 10 6 - 10" 4 Z o 
e.g. |Schneider et al.|2006| [Omukai et al.|2005| l cooling through 
line emission from heavier atoms (C,0) and molecules (H2O, 
CO, O2) starts to dominate. Thereby, the second generation of 
stars with (likely) lower masses and "Salpeter-like" initial mass 
function start to form (Population II; PopII) (see e.g. Bromm 



et al.|2001j|Schneider et al.|2002[ [2003] |Bromm & LoebT2 003 ; 



Schneider|2006| l. The transition from dominant PopIII to PopII 
star formation could already happen at early times (e.g. z » 
7), since pair-instability and core-collapse supernova explosion 
from PopIII stars can effectively enrich their environments with 
metals (Schnei der et al.|2002 Scannapieco et al.|2003| Bromm 



|& Loeb|2003HTornatore et al.|2007| >. 

Direct observations of this early period of the universe are 
challenging: halo stars with extremely low metallicities have 
been detected in our galaxy (Christlieb et al. 2002) 1, but the 
observation of a true PopIII star with zero-metallicity is still 
pending. The upcoming satellite experiment James Webb Space 



Send offprint requests to: M. Raue, e-mail: 
martin . raue@mpi -hd . mp g . de 
1 See O'Shea et al. 1 2008} for a discussion on naming conventions. 



2 For a more complete discussion on formation of the first stars 
see e.g. the 2008 updated version of |Ciardi & Ferrara|[2005| astro- 
ph/0409018 
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Telescope (JWST^] expected to be launched in 2013, with high 
sensitivity in the 1-10/mi near-infrared (NIR) band is aiming 
to detect the redshifted ultraviolet (UV) to optical (O) emission 
from source at high redshifts z > 10. 

Other constraints on the PopIII stars can be derived from in- 
tegrated properties like e.g. the number of baryons bound in stars 
or the number of ionizing photons produced. If the contribution 
from other sources to these integrated properties are reasonably 
well known, the contribution from PopIII stars can be derived. 
This can then be compared with model calculations for different 
PopIII scenarios. Tumlinson (2006) simulated the formation of 
PopIII stars using galactic chemical evolution models and com- 
pared the model output with the present day metallicity distribu- 
tion function (MDF) of the Galaxy. They found that, while not 
yet formally conclusive, the MDF could best be described by a 
PopIII initial mass functions which includes lower mass stars. 
Naga mine et al.| (2006 1 used several integrated properties in- 



cluding the extragalactic background light (EBL; see next para- 
graph) to derive constraints on the cosmological star formation 
history of PopII stars. 

In the optical to near-infrared (O-NIR) wavelength regime of 
the diffuse meta-galactic photon field (extragalactic background 
light; EBLQ stars are the main contributors to the EBL density. 
Luminous PopIII stars can leave a distinct signature in the EBL 
density ( |Bond et al.|1986] >. In particular, their contribution may 
exceed significantly the EBL density inferred from low redshift 
(z < 5) sources. Direct measurements of the EBL are difficult 
due to dominant foregrounds in our planetary system (zodiacal 
light) and the Galaxy (Hauser ~et al.| 1998] >. Nevertheless, the dis- 
covery of such a NIR background excess (NIRBE) with high 
significance has been claimed by Matsumoto et al. ( 2005|l, while 
other data showed a marginal excess (see Hauser & Dwek 2001 
for a review). The nature of this excess is still under debate. 
Dwek et al. ( |2005a[ ) find that it is likely a foreground artifact 



from zodiacal light and not of extragalactic origin. |^| A possi- 
ble PopIII origin of the NIRBE has been investigated by many 
autho rs (|Santos et al.||2002| [S alvaterra & Ferrara 2003; Dwek 
et~aT1|2005al |Madau & Silk| 20051 ISalvaterra & Ferrara(|20061 
Fernandez & Komatsu||2006| >. While |Dwek et al.| ( |2005a| > and 
Madau & Silk| ( |2005) l argue that the number of stars required 



to produce such an excess would overproduce todays metallic- 
ity and would lead to a too high number of baryons in stars, 



Fernandez & Komatsu (2006) (FK06) find that, if accounting 



for the final stage of the first stars in more detail, a PopIII origin 
of the NIRBE seems possible. 

An indirect way of deriving constraints on the EBL comes 
from the measurement of very high energy (VHE) y-ray spec- 
tra from distant sources (Ste cker et al.|1992] >. VHE y-rays inter- 
act with low energy photons from the EBL via pair-production 
( |Nikishov|19"62l|Gould & Schreder|1967] l. The cross-section of 
the pair-production is strongly peaked, so this process leaves an 
energy dependent attenuation signature in the measured VHE 
spectra. With assumption about the source physics, upper lim- 
its on the EBL density can b e derived (e.g. |Dwek & Krennrich 
2005 |Aharonian et al.|2006| |Mazin & Raue|2007| l. |Dwek et al 



( 2005b I considered the effect of a high NIRBE on the spectra 



3 http://ngst.gsfc.nasa.gov 

4 We will use the term extragalactic background light (EBL) to de- 
note the diffuse meta-galactic photon field in the UV to IR wavelength 
regime. 

3 In addition, 



Mattila 



( 2006 1 argued that the claimed discontinuity in 



the EBL at UV-0 wavelengths, which has been interpreted as a signa- 
ture for the first stars, is also an artifact of foreground subtraction. 



of distant blazars and concluded that such a high density as re- 
ported by Matsumoto et al. ( 2005 i seemed unlikely. Recently, 
strong limits on the EBL density in the NIR have been derived 
(e.g. Aharonian et al. 2006 ; Mazin & Raue 2007 ), which exclude 
the claimed NIRBE with high significance, and are only a factor 
~ 2 above the lower limits derived from source counts (Madau 
|& Pozzetti|2000] i. 

In this paper these limits on the diffuse EBL density are used 
to derive constraints on the properties of the PopIII/LM PopII 
stars. Results from a detailed model calculation of the EBL for 
different PopIII/LM PopII star scenarios are compared with re- 
cent limits on the EBL density. Our model accounts for the time 
evolution of the emissivity of a stellar population, which, for the 
case of low mass stars with long lifetimes, has profound impli- 
cations for the resulting EBL. 

The paper is organized as follows: In Sec. [2] the model cal- 
culations for the EBL density from PopIII/LM PopII stars are 
described. In Sec.[3]the resulting EBL density for different sets 
of PopIII/LM PopII star parameters is calculated and compared 
with recent limits. Limits on cosmological star formation rate 
(SFR) are derived and the detectability of a cut-off in high en- 
ergy spectra resulting from EBL attenuation are discussed. In 
Sec.|4]the derived limits are compared with previous results and 
the consequences for the PopIII/LM PopII star properties are dis- 
cussed. We summarize our results in Sec. [5] 

Throughout this paper flat Friedman cosmology is adopted 
with Qo = 0.25 , Q.\ = 0.75 and a Hubble constant of Hq = 
70 km s _I Mpc~'. 

2. EBL Model 

Previous models for the EBL density produced by PopIII/LM 
PopII stars often focused on very massive stars with M ~ 
IOOO Mq dSantos et al.|2002l|Dwek et al.|2005al|Madau & Silk| 
|2005[ ). For such very massive stars, simplified assumptions can 
be made in the calculations, e.g., the stellar emission follows 
a black body spectrum, the total emission is dominated by pho- 
tons reprocessed in the surrounding HII region and the lifetime is 
short enough to neglect the luminosity evolution. Here we want 
to explore more realistic scenarios for PopIII stars, which likely 
form with a wider range of masses. We therefore focus on a time 
dependent modeling of the evolving emissivity. 

We follow the equations presented in Kneiske et al. ( |2002 1 to 
calculate the EBL. The two basic ingredients for the calculation 
are (i) the cosmic star formation rate p* and (ii) the luminosity 
L v (t) of a specific stellar population of age t as e.g. derived from 
stellar population syntheses models. The co-moving emissivity 
(luminosity density) at redshift z is obtained from the convolu- 
tion 



£v(z) 



LMz)-t'{z'))p*{z) 



dt' 

df 



dz 



(1) 



where the cosmic star formation rate p*(z') is assumed to be- 
gin at some finite epoch z m = z(t m ) and t(z)/t'(z') is the cosmic 
time corresponding to a redshift z/z'- For given evolution of the 
emissivity a second integration over redshift yields the energy 
density, or, after multiplication with c/4n, the co-moving power 
spectrum of the EBL 



P V (Z) = VI Y (Z) = V 



4n 



dt' 
dz' 



dz' 



(2) 



with V = v(l + z')/(l + z). Further details on the calculation can 
be found in Kneisk e~et al.| (|2002 ). 
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Fig. 1. Normalized stellar initial mass functions (IMFs) for the 
IMF test cases A to E from Tumlinson (2006). All Tumlinson 
IMFs follow a lognormal distribution with peaks in the range be- 
tween 10 and 80 M Q . The values in brackets in the legend denote 
the peak position and the width of the lognormal distribution. 
For comparison a power law IMF with Salpeter slope a — 1.35 
is also shown. 



Photons from the star ionize the dense gas cloud surround- 
ing the formation site. This HII region will re-emit photons in 
emission lines, free-free and free-bound, and two-photon emis- 
sion (nebula emission). For the hot and massive stars the nebula 
emission can become dominant (e.g. Schaerer 2002). Therefore, 
two contributions to the specific luminosity are considered L v (t): 



L v (t) = Lf 



s (/) + L™ ba \t) , 



(3) 



with Z4' als (f) the stellar emission and L" ebula (f) the nebula emis- 
sion. Since we are interested in an averaged and integrated prop- 
erty, the EBL, we keep the model as simple as possible only 
considering dominant contribution to the overall emission. 

2.1. Stellar emission 

Massive PopIII stars can effectively enrich their environment 
with metals via supernova explosions (see Sec.[T]i. It is therefore 
possible, that already at early times the total stellar emission is 
dominated by the emission from PopII stars with low metallicity. 
To probe the effect of these two scenarios, we will use two stellar 
populations with different metallicities for the EBL calculation: 
(1) stars with zero metallicity corresponding to dominant PopIII 
star emission and (2) stars with low metallicity corresponding to 
dominant PopII star emission. In reality, of course, the transition 
between PopIII and PopII stars is likely an extended and patchy 
process, with parallel formation of stars in zero and metal en- 
riched environments (e.g. |Tornatore et al.|2007) >. Here we focus 
on just the two extreme scenarios to identify possible differences 
and observational signatures. 

2.1.1. Primordial/Zero metallicity (ZM) 

Emissivities for primordial/zero metallicity (ZM) stars are 
based on the stella r models from |Tumlinson et aL] ( |2003| ) and 
Tumlinson (2006). Tumlinson et al. ( 2003) ) solved the coupled 
stellar equations for ZM stars with a variant of the relaxation 
method taking into account proton-proton, CNO and He burning 
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Fig. 2. Spectra of stellar populations as function of age t. Top: 
Primordial/Zero metallicity TumA from Tumlinson ( 2006[ ); t 
ranges from to 1 Gyr (top to bottom). Bot tom: Low metallicity 
Z = 10~ 4 fromlBruzual & Chariot (2003); t ranges from to 
13 Gyr (top to bottom). 



for energy production. The stellar atmospheres have been sim- 
ulated with the TLUSTY code. Stellar population spectra (SPS) 
for a synthetic stellar cluster have been constructed from evolu- 
tionary tracks. Tumlinson (2006 ) connected present day obser- 
vational data, like the galactic metal distribution function, with 
properties of the PopIII stars via galaxy chemical evolution mod- 
els. The author tested a wide range of different initial mass func- 
tions (IMF; all lognormal, see e.g. Larson 1973 1 and constructed 
five IMF test cases labeled A to E (TumA-E), which are shown 
in Fig. [T] All IMFs test-cases are top heavy with peaks in the 
range between 10 and 80 M Q and different widths. We will adopt 
the TumA case as our fiducial model and discuss the differences 
in the EBL resulting from the different IMFs in Sec. |3.1| The 
SPSs for the case TumA and different ages are shown in Fig. [2] 
top panel. Nuclear burning in ZM stars occurs at higher temper- 
atures than in PopII stars (Tumlinson et al. 2003 ), which results 
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in a large number of H and He ionizing photons in the UV. High 
mass stars have a short lifetime, so the overall emission of a high 
mass stellar population is significantly reduced at stellar popula- 
tion ages > 0.1 Gyrs. 



2.1.2. Low metallicity (LM) 

For stellar populations with low metallicity (LM) the well es- 
tablished Isochrone Synthesis Spectral Evolutionary Code from 
IBruzual & Chariot (2003 ) is utilized. This code has also been 



used for the PopII star component in the EBL model of Kneiske 



et al. 



( |2002fr . Here we choose a metallicity of Z = 10 and a 
Salpeter IMF. In contrast to the lognormal IMFs of the ZM stars 
described in the previous section, this will lead to a large num- 
ber of stars with masses < 10 M Q (Fig.[T]i. The SPSs for different 
ages are shown in Fig.|2]bottom panel. Due to the longer lifetime 
of the low mass stars, the overall emission time of the LM pop- 
ulation is significantly longer than in the case of ZM stars with 
high masses. 

The LM second generation stars already form in (low) dusty 
environments. The dust formation at such high redshifts is 
not known and strongly depends on the environment. Dust re- 
emission would occur in the MIR to FIR and redshifts into the 
microwave regime, which is completely dominated by the cos- 
mic microwave background (CMB). For the wavelength range 
considered here (O-NIR) dust absorption is negligible and we 
will not consider it further. 



2.2. Nebula and line emission 

The PopIII stars are embedded in dense gas clouds, which ab- 
sorb and re-process photons from the star, creating an H II region 
(nebula). They burn at higher temperatures than their successors 
( Tumlinson et al. 2003[ ), producing copious amounts of H and He 
ionizing photons. The nebula absorption and emission can there- 
fore severely alter the SPS spectra ( Santos et al.|2002| |Schaerer 
[20021 [Fernandez & Komatsu|2006| ). 

In our calculations, we assume total absorption of hydrogen- 
ionizing photons in the nebula. As discussed in |Fernandez & 
|Komatsu] ( |2006| l, the continuum luminosity does not depend on 
the number density of electrons or protons in the nebula, since, in 
case of a Stromgren sphere, a higher number densities would re- 
sult in a higher re-combination rate and consequently a smaller 
emission region. This approximation is correct as long as the 
bulk of the emission comes from the nebula around the star in 
the host halo. The effect of escaping UV photons into the inter- 
galactic medium (IGM) on the SPS spectrum has been calculated 
by Santos et al. ( 2002 ). The absorption of escaping ionizing pho- 
tons would then take place in the IGM leading to a similar spec- 
trum in the O-NIR wavelength range as in the case of absorption 
in the nebula0 

To keep the calculation simple, we will only consider the 
continuum and line emission of hydrogen, following the calcu- 
lations of Fernandez & Komatsu (2006) (FK06 in the follow- 
ing). The contributions from free-free and free-bound continuum 
emission, 2-photon emission, and the Lyman-o' line emission are 
taken into account. Below, the equations governing the nebula 
emission are briefly summarized. For details we refer the reader 
to FK06. 



6 At larger wavelengths, the spectra are different due to the different 
strength in the free-free continuum emission. 
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Fig. 3. Stellar and nebula emission for stellar population of 
age 10 3 (upper panel) and 2 x 10 7 years (lower panel) (stel- 
lar emission: ZM, TumA). The emission components shown 
are stellar emission (green solid), Lyman-o- line (red dashed 
double-dotted), 2-photon (purple dashed), free-bound (light-blue 
dotted), and free-free continuum emission (dark-blue dashed- 
dotted). The thick solid line is the total emission. 



The luminosity of the free-free and free-bound continuum 
emission is given by 



j^cont _ 



n e n p a 



(4) 



where Q H = j°° L s v ' ar v~ l dv is the production rate of hydrogen 
ionizing photons, n e and n p number density of electrons and pro- 
tons, respectively, ag ~ 2.06 ■ 10 T, 



-1/2 3 

cm 



s 1 is the case- 



B recombination coefficient for hydrogen at a temperature of 
T g = 20.000 K (see |Spitzer|1978| . 

The total emissivity for free-free and free-bound emission 



taken from Dopita & Sutherland ( 2003 ) is 



e v — Ann, 



n p r c 



exp(-hv/kTg) 

T l/2 



(5) 
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The coefficient y c is a constant Fk '- 
free-free and free-bound emission 



5.44-10 39 times a term for 



Fkiigff) + £ 



(gfb)) 



(6) 



with x n = Ry/(kT g n 2 ) and (g/f) ~ 1.1 and {gfb) ~ 1-05 are the 
gaunt factors for free-free and free-bound emission, respectively 
(FK06). Note that for free-bound emission the sum is over all 
bound states with energies below the frequency in question and 
the infinite sum over n. 

The luminosity for the line contribution is given by 



Lf e = Y J hv i cf> i (v-v i )f i Q H 



(7) 



where </> ; (v - v,) is the line profile which is assumed to be a delta 
function. The fraction of ionizing photons which are converted 
to a line i is fi. To keep the model simple only the Lyman-a 
line emission will be included, since the effect of other series is 
very small ( |Santos et al.||2002| |Dwek et al.|[2005a| l. A fraction 
°f fiya - 0.66 is used (see Spitzer|1978 for 20.000K), because 
2/3 of the time a n=2 — > 1 transition creates a Lyman-o- photon 
and 1/3 of the time a continuum photon via 2-photon decay is 
emitted. 

The luminosity of the two-photon emission is calculated by 



I%= — (l-f Lyci )P(v/v lya )Q H 

v Ly a 



2hv 



(8) 



where P(y)dy is the normalized probability per 2-photon decay 
of getting one photon in the range dy = dvlvi Ja . An analytical 
expression for P(y) has been derived from a polynomial fit to the 
data given in Table 4 of Bro wn & Mathews] (|1970| l, which also 
fits well the data for y > 0.5 (e.g.|Gaskell|1980)i: 



P(y) 



1 .307 -2.627 '(y 
-51.69(y 



0.5) 2 +2.563(y-0.5) 4 



(9) 



0.5) 6 . 



The different components of the nebula emission together 
with the stellar and total emission for a stellar population of age 
10 3 and 2 X 10 7 years are shown in Fig. [3] (ZM stellar popu- 
lation model, TumA). For young stellar populations the neb- 
ula emission dominates the emission spectrum at wavelengths 
A > 1200 A. Older stellar populations produce less ionizing pho- 
tons, so the nebula emission gradually decreases until it is negli- 
gible. 

2.3. Star formation rate 

The star formation rate (SFR) at high redshift (primordial SFR; 
PSFR) is largely unknown. Direct measurements and limits on 
the PSFR from deep surveys of UV-bright galaxy with the 



et al. 



2006 



2007 



dropout technique (Bouwens & Illingworth 2006 ) |Richard et al 
IMannucci et al.|2007HBouwens et al.|2007||2008HRichard 



2008| l, Ly-q emitters pCashikawa et al.||2006| |Ota et al 



2008 , and gamma-ray burs ts (|Firmani et aL 2004 Le & Dermer 



|Guetta & Piran|2007|[Yuksel et al.|2008| ) exist, but the in 



dividual measurements do not agree. While some measurements 



show a strong decline in the SFR at z > 5 (e.g. Mannucci et al 
2007 Bouwens et al. 200811, others detec t a constant or only 



slowly dec lining SFR (e.g. [Richard et al.|[2006l [20081 [Yiiksel 
et al.|2008) >. The different methods most likely sample different 
populations of PSFR regions and are not directly comparable. 
The individual measured PSFRs therefore can be considered as 
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Fig. 4. Sketch of the parameterization for the star formation rate 
(SFR) as broken power-law (Eq. [T0j >. A set of SFR parameter 
will be denote as SFR(z pea k, a, ft), with the normalization set to 
Po = 1- 

a lower limit to the total PSFR d Yiiksel et al. 12008b. Furthermore, 



PopIII stars form in small (proto-galactic) halos (e.g.|Greif et al. 



2008 ), so the galaxy PSFR, as derived from UV bright galaxies, 
may not be representative for the PSFR. 

Further constraints on the PSFR come from constraints on 
the reionziation history of the universe. The detection of absorp- 
tion features in the spectra from distant quasars (Gunn-Peterson 
trough; Gunn & Peterson|1965| l can provide direct evidence on 
the ionization state of the early universe. Such observations in- 
dicate that the universe was largely ionized at z ~ 6 (Fan et al. 
2006 ). Recent measurements of the fluctuation power spectrum 
and polarization of the CMB from the WMAP 5 year data also 
point toward an early and extended reionization, with the uni- 
verse mostly ionized at redshift z re ion ~ 1 1 ± 1.4 (Dunkley et al. 
2008 ). While the sources of the reionization are not well known, 
PopIII/LM PopII stars are considered to be the natural candi- 
dates HBarkana & Loeb|2001| i. 

For our model, we adopt a simple broken power law to de- 
scribe the primordial SFR: 



P*(z) = Po 



z + l 



(10) 



with q - a for z < z pea k and q = ft for z > z pea k (sketched in 
Fig. [4]). We only consider star formation in the redshift range 



from z start = 35 to z 



stop 



5, i.e. p*(z > z start ) 



and 



p*(z < Zstop) = 0. In the following a set of SFR parameters will 
be denoted as SFR(z pea k, a, ft). 

2.4. Resulting EBL 

The EBL density for the ZM and the LM model for redshifts 
z = 8, 4, are shown in Fig. |5](SFR(10, 10, -2)). For high red- 
shifts, when stars are still being actively formed, the EBL density 
from ZM stars is factor 2-3 higher than in the LM case. This in- 
verts for low redshifts: at z = the EBL density from LM stars 
is factor ~2 higher than in the case of ZM stars. The high mass 
ZM stars have all ended their lifetime and no new photons are 
produced, while the low mass LM stars are still burning. The 
contribution from low mass stars also blurs the sharp EBL den- 



ZM 



LM 
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Fig. 5. Co-moving EBL energy density for ZM stars (top) and LM stars (bottom) for redshift z = 8 (left), z = 4 (middle), 
(right) and SFR(10, 10, -2). The contribution from the nebula (red dashed) and the stellar (blue) components are shown. 



sity drop-off at the ionization edge in the UV, resulting in an 
overall broader shape and more low wavelength photons than 
for ZM stars. Consequently, no sharp drop in the EBL at < 1 fim 
is expected for the case of LM stars. For the hot and massive 
ZM stars the nebula emission makes the dominant contribution 
to the EBL at the longer wavelengths. For the LM model the 
contribution from the nebula is negligible. 

These results can be compared with the results from FK06, 
which only considered stars with mass > 5 MoQ they find that 
the nebula and line emission are always the dominant contribu- 
tors to the EBL density, even in the case of low metallicity stars 
with Salpeter IMF (FK06, Tab. 1). When accounting for the low 
mass range of the IMF (< 5 M ), the emission from old stars 
dominate the EBL and the nebula contribution is negligible. 

3. Constraining the early star formation 

As discussed in Sec. [TJ the EBL can be used as probe for the 
properties of the first stars. While several authors discussed a 



possible signature from the PopIII stars in the NIRBE (Santos 
eTaTI [20021 |Salvaterra & Ferrara| [20031 |Dwek et~ai"l|2005al 
Mada u & Silkl|20051 |Salvaterra & Ferrara||20061 |Fernandez & 
Komatsu 2006 ), we want to focus on the recent limit on the EBL 
density ( |Aharonian et al.|2006[ |Mazin & Raue|2007| l and what 
constraints on the properties of the PopIII/LM PopII stars can be 
derived from these limits. Two properties will be mainly investi- 
gated: (1) the SFR and (2) the influence of the metallicity. 

The normalization of the SFR p ( > (at z = z pea k) directly trans- 
lates into the overall normalization of the EBL (see Eq.[TJand[2]i. 
We will set po = 1 throughout our calculations and later use the 
scaling relation to derive constraints on po- The SFR parameters 



7 More precisely, stars with lifetimes short enough so that the approx- 
imation used by FK06 (no time evolution of emissivity, only averaged 
quantities) is correct. 



Parameter 


Values 


£peak 


7, 8, 9, 10, 11, 12, 13, 14 


a 


4, 7, 10 




-2,-1,0 



Table 1. Parameter-values of the broken-power-law SFR 
(Eq.m. 



tested are given in Tab. [TJ The choice of z pea k is motivated by 
the limits on the redshift of reionization derived from WMAP 
data, which gives a 3 <x lower limit for a sudden reionization of 
Zreion > 6.7 with a best fit of Zreion ~ 11 + 1.4 and evidence for 
an extended process (Dunkl ey et al.|2008] l. A steep slope a = 10 
corresponds to a sudden stop of the PopIII star formation and 
negligible contributions at low z while flatter slopes (a = 4) give 
some contribution of the PopIII SFR even at low z as, e.g., pre- 
dicted by |Tornatore et aL (|2007) from numerical simulations. 
/3 — corresponds to an extreme case of a constant SFR (up to 
Zpeak) (FK06), while f3 = -2 corresponds to a more standard case 
(e.g. |Dweketal.|2005a| >. 

Figure [6] shows the impact of the different SFR parameters 
on the resulting EBL energy density. Shifting the peak position 
Zpeak results in a shift in the wavelength of the maximum of the 
EBL density, and a lower EBL density for higher z (first row). 
Changing the slope a has a similar effect as changing the peak 
position, i.e. shifting more of the star formation to smaller/higher 
z values (second row). The slope ft changes the EBL density at 
wavelength greater than the peak in the EBL density (third row). 

The EBL density from LM stars (3rd row) has a broader 
peak extending to lower wavelengths and a slightly higher over- 
all EBL density than the EBL density from ZM stars (2nd row), 
mainly due to the longer burning times of the low mass stars. If 
the first stars have a significant contribution to the overall EBL 
in excess of the contribution from PopII stars, the detection of 
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Fig. 6. Co-moving PopIII SFR (left) and resulting EBL energy density (ZM: middle; LM: right). In each row one parameter of the 
SFR parameterization is varied. As reference for the low redshift SFR the data collection from Hopkins & Beacom (2006) is given 
(grey markers). Top: z pea k varies - SFR(x, 10, -2). Middle: a varies - SFR(10, x, -2). Bottom: B varies - SFR(10, 10, x). 



a break in the EBL from UV to NIR could be interpreted as an 
indication for massive ZM stars with short lifetimes (e.g. [Santos] 
et al. 2002 ). Furthermore, the EBL from LM stars shows a differ- 
ent redshift dependency (Fig.|5]l, which in principle will result in 
a different absorption signature for high redshift VHE sources. 
Detecting such a signature will be a difficult task: high precision 
measurements of VHE spectra for high redshift are needed and 
more important, the PopIII contribution to the EBL density has 
to be in excess of the contribution from second generation of 
stars (see Sec.|3.2|for a detailed discussion). 



3. 1 . Initial mass function 



The zero metallicity SPS from Tumlinson (2006) are calculated 
for several different IMF cases shown in Fig. [1J Figure [7] shows 
the EBL density resulting for the different IMFs TumA-E using 
the same SFR(10, 10, -2). For the wavelength range of inter- 
est (UV to NIR; < 4 /mi), the differences are very small (or- 
der < 10%). To differentiate between such subtle differences the 
EBL density has to be resolved on the same level. As discussed 
below (Sec. |3.2[ ), current EBL limits constrain the EBL density 
in this wavelength range in the order of < 50% at best, so for 
further calculation we will only consider the TumA IMF case 
(somewhat average EBL) for the ZM stars. 



10 



Iff 1 : 



Iff 1 




X (p.m) 



Fig. 7. Resulting EBL energy density for different IMFs 
(Tumlinson A, B, C, D, and E) for SFR(10, 10, -2). Given the 
similarities between the EBL densities, with only small differ- 
ences at longer wavelengths, we will use the TumA case for fur- 
ther calculations. 
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Fig. 8. EBL density at 2 fim (z - 0) versus the peak of the SFR 
Zpeak f° r different stellar models and SFR parameters. Left: ZM. 
Right: LM. Other SFR parameters: a — 10 (solid lines), a = 
4 (dashed lines); ft = (triangles), /3 = — 1 (circles), (3 — -2 
(squares). Results for a = 7 lie between a — 10 and a — 4 and 
are omitted to improve the readability of the plot. 




Fig. 9. Upper limits on the SFR (po) derived from the EBL 
contribution at 2/im assuming a maximum EBL contribution 
of 5nWm _2 sr'. Models: ZM (black lines and markers), LM 
(red/grey lines and markers). SFR parameters: a = 10 (solid 
lines), a = 4 (dashed lines); /3 — (triangles), j3 — -2 (squares). 



3.2. Constraints from HE/VHE observations 
VHE spectra 

The contribution from a redshifted early stellar population to 
the EBL density is expected to peak in the wavelength range 

( 2006[ ) derived up- 



around 1-2 jjm. Recently, |Aharonian et al.| 
per limits of v/ v < (14 + 4)nWirT / sr" 1 on t ' 
this wavelength range, utilizing the hard VHE spectrum of the 
blazar 1ES 1101-232 (located at a redshift z = 0.186) and cer- 



le EBL density in 



tain assumptions about the sour ce physics (see also|Katarzyriski 
let al.|20061|Stecker et al.|2007l|Aharonian et al.|2008| for some 
caveats). This limit has been confirmed by the detection of a 
second source of VF1E y-rays with similar properties 1ES 0347- 
121 (z = 0.188) ( |Aharonian et al.||2007] l and was shown to 
be largely independent of the exact shape of the EBL density 



flMazin & Raue|2007| l. pvladau & Pozzetti| (|2000) derive a strict 
lower limit on the EBL density of > (9.0 - 9.7^'°) nW nr 2 sr~' 
from deep source counts (not corrected for completeness) in the 
wavelength range 1 - 2 fim. [Totani et ah\ ([2001 ) estimated the 
EBL density contribution from resolved galaxies to be (10.1 - 
12.8)nWm- 2 sr 1 at 1.25 nm and (7.8 - 10.2)nWirT 2 sr 1 at 
2.2/im, accounting for missed galaxies due to selection effect. 

In this paper we will adopt a limit on the EBL contribution 
from the PopIII/LM PopII stars of ~ 5 nW irT 2 sr -1 at 1 - 2 ji/mj^] 
Comparing this maximum contribution with the EBL density 
value at 2 fim Febl(2 yum), as calculated from our model for a 
specific set of parameters, an upper limit on the normalization of 
the SFR po can be derived: 



p <5nWm sr /F EBL (2fim) 



(11) 



whereas Febl(2 //m) is calculated using po = 1 M G Mpc" 3 yr~'. 

The resulting EBL densities at 2/im for different SFR 
parameter-sets and metallicities are shown in Fig. [8] The 
EBL contribution for the ZM model range from 1.5 to 
8.5nWm -2 sr 1 , for the LM model from 3 to 15 nWrrT 2 sr -1 . 
Converting these EBL contributions to limits on the normaliza- 
tion of the SFR, limits from 0.3 to 3 M Q Mpc~ 3 yr~' are derived 
(Fig. |9j. LM models result in an overall factor ~2 stronger con- 
straints due to the higher EBL contribution. The different SFR 
parameters result in a similar spread in the limit of factor 2-3. 
The limits on the SFR directly scale with the limit on the EBL 
density contribution at 2/im: if e.g. the EBL limit is lowered 
from 5 to 1 nW irT 2 s~' the corresponding SFR limit is also low- 
ered by a factor of 5. 



Cut-off in high redshift sources 

For distant source (z > 1) of HE/VHE y-rays, the attenuation 
from the EBL results in a strong suppression (cut-off) at ener- 
gies > 10 - 100 GeV. The cut-off energy shifts with redshift, 
which in principle could provide a measurement of the EBL den- 
sity ( Fazio & Stecker|1970 1. The cut-off is expected at energies 
< 100 GeV, which will be probed by the Large Area Telescope 
(LAT) instrument on the recently launched FERMI satellite (en- 
ergy range: ~20MeV - 300 GeV). 

The main challenge in determining such a cut-off is to dis- 
tinguish it from source intrinsic effects like e.g. insufficient ac- 
celeration of particles to high energies or internal absorption, 
which could also lead to a cut-off at high energies. Source intrin- 
sic absorption can even lead to a redshift dependent absorption 
signature, mimicking the z dependent EBL attenuation (Reimer 
|2007| ). Furthermore, energies ~ 100 GeV are at the outer edge of 



the LAT sensitivity band, with the sensitivity decrease by more 
than one order of magnitude in comparison with the best sensi- 
tivity around 1 Ge"V|_J 

Detecting an absorption signature from the PopIII/LM PopII 
stars is even more challenging: in addition to source intrinsic 
effects, the absorption from the PopII EBL acts as foreground, 
from which the absorption from the PopIII EBL has to be dif- 
ferentiated. Figure 10 shows the optical depth derived for the 
ZM PopIII model in comparison to the optical depth from the 
low PopII model from Kneiske et al. (2002 2004| for a source 



8 This corresponds to the best guess, not the extreme limit, i.e. using 
the outer error range, which is of order ~ 9 nW m 2 sr~' . 

9 http://www-glast.slac.stanford.edu/software/IS/ 
glast_lat_performance.htm 
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Fig. 10. Optical depth for VHE y-rays derived for the ZM model with SFR(7 - 14, 10, -2) (black solid lines) versus the optical depth 
derived for a PopII EBL model (low model from Kneiske et al. 2002, 2004; red dashed line) for a source at redshift z — 1 (left) and 
z = 2 (right). 



at redshift z = 1 and z - 2[^]For redshift z — 1, the optical 
depth resulting from the PopII EBL dominates over the PopIII 
contribution (with po = 1, which is of order of the limit derived 
in the previous section). For higher redshifts the situation is dif- 
ferent: the PopII SFR rises steeply with redshift with a plateau 
or peak expected in the redshift range z = 1 - 2, so the main 
PopII EBL contribution is building up over the redshift range 
Z = — 2. The PopIII will not add new photons to the EBL in 
this redshift range, so the co-moving EBL density contribution 
from PopIII stars is constant and just scales with the cosmolog- 
ical expansion/contraction. Consequently, at redshift z > 2 the 
optical depth resulting from PopIII EBL can dominate the total 
attenuation at energies >30GeV (Fig. 10 right). For a source at 
redshift z ~ 2 the optical depth resulting from the PopII EBL is 
> 1, so to distinguish between PopIII and PopII contribution, one 
would have to differentiate between optical depths of e.g. 1 and 
10. Note that the PopII EBL model could underestimate the EBL 
density in the UV-O range, since the contribution from AGN is 
not included in the calculation. Thus the background resulting 
from other contributors to the EBL could be even higher. 

In case of the LM scenarios the resulting optical depth is 
higher than for the ZM models, due to the higher EBL density 
at A < 0.5/mi (Fig. [6j, but this will not change the fundamental 
challenge of determining the exact shape of a steep cut-off from 
a low statistic measurement. 



4. Discussion 

Figure[TT]shows the limit on the PopIII SFR derived in this paper 
in comparison with other measurements and limits on the SFR. 
As argued in Sec. |2.3| the direct measurements (from e.g. source 
counts) are not directly comparable with the limits derived from 
an integrated property like the EBL, since they most likely sam- 
ple different contribution to the total SFR. In this respect, they 
have to be considered lower limits on the total SFR. The direct 




Fig. 11. Limits on the co-moving SFR of PopIII stars derived in 
this paper in comparison to other limits and measurements of the 
star formation rate. The purple striped region indicates the range 
given in FK06 (see text for details). Blue open markers are from 
the data collection from Yii ksel et al.| ( |2008| l, while blue filled 
markers at z = 4.5 and z — 6 are the data points derived in the 
same paper from GRB measurements. As reference for the low 
redshift SFR the data collection from Hopkins & Beacom (2006) 
is given (grey markers). 



measurements at z ~ 7 lie two to three orders of magnitude be 
low our limits, while rec ent determination of the SFR at z = 3-6 
give 0.1 Mq Mpc~ 3 yr _I (Yiiksel et al 



2008) 



Henry et al. (2008) 



reported a candidate for a Lyman break galaxy at z = 9. If con- 
firmed, this could imply that the SFR at z ~ 9 does not evolve 



10 For details on the calculation of the optical depth sec e.g. 
R aue|2007 



Mazin & 



strongly and is roughly on the same level as at z = 3 ( |Henry etaL 
2008 ), which would bring our limits in range (factor 5-10) of the 
direct measurements. 
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Other results on the SFR of PopIII stars come from numer- 
ical simulations: Brom m~& Loeb| ( |2002) > derive a peak SFR 
of 0.3-0.6 M Mpc~ J yr -1 depending on the dominant cooling 
mechanism. Tornatore et al. (2007) find in their simulation low 
SFRs with peaks of 10" 3 andlO^ M G Mpc 3 yr 1 at z * 6 for 
PopIII and PopII stars, respectively. 

Valuable constraints on the properties of PopIII star have also 



been derived from reionization studies. |Choudhury & Ferrara 
(2006) used a self-consistent reionization model together with 



experimental results to derive a best fit model of the reion 
ization history. Their best fit model, which predicts negligible 
source counts at z > 10, has a peak SFR for the PopIII stars 
of order 10~ 3 M Q Mpc~ 3 yr (z ~ 12) and is always dominated 
by LM PopII stars (PopII SFR rising to 0.1 M Q Mpc 3 yr _I at 
z — 7, smoothly connecting with the low redshift PopII SFR). 
However, if some of the high redshift source candidates from 
deep surveys (see Sec. 2.3 i would turn out to be valid z ~ 10 
sources, requiring a higher SFR at z = 10, the model would 
be strongly constrained by the limits on the Ly a optical depth, 
possibly allowing to pin down how reionization occurred. Greif 
& Bromml ( |2006| ) investigated the reionization history with a 
semi-analytic model utilizing two distinct ZM PopIII star popu- 
lation^] They investigated the effect of several feedback mech- 
anisms (Lymann-Werner, photoheating, chemical feedback) on 
the SFR of PopIII stars and found the peak SFR of PopIII in 
the order of 10" 4 - 10 3 M Q Mpc 3 yr -1 . Again, the overall SFR 
in the redshift range z = 10 — 15 is dominated by ( LM) PopII 
stars (~ 10" 2 - 10" 1 M G Mpc 3 yr 1 at z = 5 - 15). |Haiman & 



Bryan ( 2006[ ) found evidence for a suppression of the SFR in 



high-redshift minihalos, utilizing the WMAP constraints on the 
optical depth to electron scattering. For zero metallicity PopIII 
stars the limits on the SFR derived in this paper are 1 to 5 orders 
of magnitude above the predictions from reionization models. 
For LM PopII the limits are still above the predictions, but only 
a factor 3 to 10. 

FK06 derived limits on the PopIII star formation rate using 
a simplified EBL model, assuming a PopIII EBL contribution of 
2 - 50nWm~ 2 sr _1 . We convert their result (FK06, Tab. 1) using 
a maximum PopIII EBL of 5 nWirT 2 sr 1 . The resulting limit 
on the PopIII SFR is 0.63 - 1.2M Q Mpc 3 yr' (Fig, [TTj purple 
striped box), which is in the range of the limit derived in this 
paper. Note that these results are not directly comparable with 
our results since (1) FK06 used a simplified model for the stellar 
emission, which does not account for the temporal evolution of 
the emission. While for ZM stars with top-heavy IMFs this is a 
valid approximation (high mass corresponds to short lifetime), 
for LM stars with Salpeter IMF the EBL contribution from low- 
mass and therefore long-lived stars changes the shape of the EBL 
and increase the overall EBL density (see Sec. 2.4 1. (2) FK06 
only consider a constant SFR similar to our case SFR(x, 10 -2). 
Considering different SFRs introduces a redshift dependency of 
the limit and a spread of factor ~ 5. 

An upper limit on the baryonic mass density bound in PopIII 
stars M, can be derived by integrating the SFR over time 



M t < AC 



lim 



f 

Jo 



(12) 



which can be compared with estimates for the total density of 
baryons locked up in stars of 6%+2% £2b (Fukugita & Peebles 
2004), with fib being the baryonic density of the universe in 



11 Classical high mass > 100 M Q PopIII stars from clouds with molec- 
ular hydrogen and lower mass > 10M o PopII.5 stars, formed in clouds 
with HD cooling. 



TumA 


TumB 


TumC 


TumD 


TumE 


Salpeter 


./■em 0.27 
./metal 0.12 


0.12 
0.15 


0.56 
0.07 


0.37 
0.11 


0.43 
0.27 


0.6 
< 0.047 



Table 2. Remnant mass fraction / lem (Eq. 13) and metal mass 
fraction / me tai for different IMFs. For / rem for the Salpeter IMF 
results from Nagamine et al. ( 2006 ) have been used. 



units of the critical density. For our limiting SFRs we derive 
Mj im *5% Q b for the LM and M^ m «10% Q. b for the ZM stars. 
Since not all baryons remain locked up in stars (they get ejected 
by supernova explosions or stellar winds) this number has to be 
corrected for baryons leaving the stars. Nagamine et al. (2006) 
calculated the fraction of amount of recycled gas to the total 
amount of gas initially converted to stars to be / rec = 0.32 for 
Salpeter IMF and solar metallicity. Adopting a remnant mass 
fraction of 0.6 of the initial mass the value for LM stars is low- 
ered to M* ~3%Qb- For the ZM stars with lifetimes ~ 10 6 yr 
it is reasonable to assume that all ZM stars created at redshifts 
z > 5 have reached their final state. We therefore estimate the 
remnant mass fraction / rem using the final state masses for ZM 
stars presented in Fig. 2 of Heger & Woosley (2 002[ ): 



/rem — 



J mm . n N(m)M rem (m)dm 
f™""* N(m) m dm 



(13) 



with N(m) being the IMF and M lem (m) the remnant masses. The 
derived fractions range from 0.12 to 0.56 depending on the IMF 
(Tab. [2]), which results in a baryon density locked up in ZM stars 
of ~ 1.2 - 5.6% Qb- The final states of ZM stars are not fully 
understood, so these numbers should be considered estimates not 
precise calculations. However, these estimates demonstrate that 
our SFR limits are below the range what could be excluded from 
the limits on the density of baryons in stars. 

Similar calculations can be performed for the metal enrich- 
ment from the PopIII/LM PopII stars. The metal mass fraction 
/metal can be calculated from Eq. [I3]by replacing M lem (m) with 
M meta i(m), the metal yield produced by a star of mass m. Here 
it is again implicitly assumed that all star have reached their fi- 
nal stage. For ZM stars this is a valid assumption, but for the 
case of LM stars this calculation results in an upper limit on the 
metal density produced by the stars. We adopt metal yields from 



Heger & Woosley (2002) (Z = 0; M > 12OM ),|Portinari et al. 



( 1998} (Z = 0.0004; 5 M G < M < 120 M Q ) and Marigo|| |200l7 T 
Z = 0.004; 0.8 M G < M < 5 M Q ). The results are summarized in 
Tab. [2 LM stars produce metal densities of < 0.2% which 
is well below todays metal density of ~ 2%Qb. ZM stars pro- 
duce metal densities of ~ 0.7 - 2.7% Qb depending on the IMF. 
Again, we stress that these are estimates, since there are large 
uncertainties in the final states of the ZM stars, especially in the 
high mass range. Still, ZM stars with very heavy IMFs (e.g. the 
Tum-E case) seem to overproduce the metal density, if a SFR on 
the level of our limits is assumed. 



5. Summary & Conclusions 

We investigate how limits on an integrated present-day observ- 
able, the EBL, can be used to constrain the parameters of the 
early stars. A detailed model for the PopIII/LM PopII star emis- 
sion from a large range of different scenarios is used to calculate 
the evolving EBL from these stars, taking into account the time 
evolution of the emissivity and the emission from reprocessed 
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ionizing photons (nebula). Recent limits on the EBL density de- 
rived by Aharonian et al. ( 2006 1 from the detection of hard VHE 
y-ray spectra from distant sources together with lower limits 
from source counts (Madau & Pozzetti|2000||Totani et aLj2 001 ) 
suggest a maximum PopHI EBL contribution of ~ 5 nWm sr 
at 1 - 2 //m. Comparing this contribution with our model calcu- 
lations, a limit on the co-moving SFR of PopIII stars of 0.3 to 
3M Mpc _3 yr _1 is derived for the redshift range 7 - 14. This 
limit depends on the redshift, on the exact shape of the SFR and 
on the assumed scenario for the early star formation: if the early 
star formation is dominated by second generation stars with low 
metallicity, the limit is factor two lower than in the case of zero 
metallicity stars. The SFR limit directly scales with the assumed 
PopIII EBL contribution, e.g., if the EBL limit is lowered by a 
factor 2 the corresponding SFR limit is also lowered by the same 
factor. 

The SFR at redshift > 5 is difficult to access via direct ob- 
servations. A few measurements and limits exist, generally fa- 
voring a lower SFR in the range of 10~ 3 to 10~ 2 M Mpc~ 3 yr -1 , 

Recent mea- 
3-6 favor 



but the spread and uncertainties are large (Fig. 1 1 
surements of the SFR in the redshift range z ■ 
a flat (or even increasing) SFR of the order of 0.1M o Mpc~ 3 
( jYiiksel et aLl|2008t |Faucher-Giguere et al.||2008) , which is in 
the range of our best limit (0.3 Mq Mpc~ J yr" 1 at z = 7 for 
LM, SFR(7, 4, 0)). Stringent constraints on the properties of 
the first stars come from reionization studies, which combine 
complex semi-analytical modeling with limits on the reioniza- 
tion history. Predictions for the peak SFR for PopIII and (LM) 

-3 , 



PopII are in the range of 10 — 10 M Mpc yr and 
1Q- 2 - lQ-'MpMpc " 



2006 Greif & Bromm 



yr ', respi 
|2006 )F1- 



, respectively ( |Choudhury & Ferrara 
While for PopIII stars the lim- 



its derived in this paper are 1 to 5 orders of magnitude above 
the SFR expected from the best fit models from these studies, 
for (LM) PopII stars the limits are close (factor 3 to 10) to these 
predictions. 

Pair-creation of VHE photons from distant sources (z > 1) 
with the low energy photons from the EBL results in a sharp 
cut-off in energy spectra > 30 GeV, which should be detected by 
the FERMI experiment. To derive constraints on the PopIII/LM 
PopII stars from the detection of such a cut-off is challenging 
since (a) the photon statistics will likely be low, (b) attenuation 
from the PopIII/LM PopII EBL competes with the attenuation 
due to the EBL from PopII stars, which is likely the dominant 
contribution to the total attenuation, and (c) the general problem 
to discriminate between source intrinsic effects and attenuation 
from the EBL. 

Constraints on the EBL can provide additional insides in the 
star formation processes of the early universe. Though the lim- 
its are not (yet) strongly constraining, they provide an indepen- 
dent probe for the star formation at redshift z > 5. With the 
current limits on the EBL in the NIR it is not possible to dis- 
tinguish between different PopIII IMFs or metallicity scenarios. 
In the future, the Cherenkov Telescope Array (CTATjwill pro- 
vide sensitive measurements in the ~20GeV to lOOTeV energy 
range, which will result in strong constraints on the EBL in a 
wide wavelength range. Together with direct detections and deep 
source counts from upcoming satellite and ground-based tele- 
scopes this will enable to resolve many of the contributors to the 
EBL and thereby tighten the limits on the PopIII/LM PopII stars 
properties derived from the EBL. 

12 A higher SFR is also possible, see Sect. [4] 

13 http://www.cta-observatory.org/ 
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